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The asteroid belt is the leftover of the original planetesimal population in the inner solar 
system. However, currently the asteroids have orbits with all possible values of eccentricities 
and inclinations compatible with long-term dynamical stability, whereas the initial planetesimal 
orbits should have been quasi-circular and almost co-planar. The total mass now contained in 
the asteroid population is a small fraction of that existing primordially. Also, asteroids with 
different chemical/mineralogical properties are not ranked in an orderly manner with mean 
heliocentric distance (orbital semi major axis) as one could expect from the existence of a radial 
gradient of the temperature in the proto-planetary disk, but they are partially mixed. These 
properties show that the asteroid belt has been severely sculpted by one or a series of processes 
during its lifetime. This paper reviews the processes that have been proposed so far, discussing 
the properties that they explain and the problems that they are confronted with. Emphasis 
is paid to the interplay between the dynamical and the collisional evolution of the asteroid 
population, which allows the use of the size distribution or of the crater densities observed in the 
asteroid belt to constrain the dynamical models. We divide the asteroid belt evolution into three 
phases. The first phase started during the lifetime of the gaseous proto-planetary disk, when 
the giant planets formed and presumably experienced large-scale migrations, and continued 
after the removal of the gas, during the build-up of the terrestrial planets. The second phase 
occurred after the removal of the gaseous proto-planetary disk and it became particularly lively 
for the asteroid belt when the giant planets suddenly changed their orbits, as a result of a mutual 
dynamical instability and the interaction with the trans-Neptunian planetesimal disk. The third 
phase covers the aftermath of the giant planet instability, until today. 


1. INTRODUCTION 

The asteroid belt helps us in reconstructing the origin 
and the evolution of the Solar System, probably better than 
the planets themselves. This is because the asteroid belt 
provides several key constraints that can be used effectively 
to guide the development, the calibration and the valida¬ 
tion of evolutionary models. Compared to other small body 
populations, such as the Kuiper belt or the Oort cloud, the 
constraints provided by the asteroid belt are probably more 
stringent, due to the fact that the number and the properties 
of the asteroids are better known, thanks to ground based 
observations, space missions and meteorite analysis. 

The structure of this review chapter is therefore as fol¬ 
lows. We start by reviewing in section 2 what the most 


important observational constraints on the asteroid belt are 
and what they suggest. Then, in section 3, we will review 
the main models proposed, from the oldest to the most re¬ 
cent ones, and from the earliest to the latest evolutionary 
phases they address. In Section 4, we will discuss several 
implications for asteroid science from our current preferred 
view of the dynamical evolution of the asteroid belt. 

The dynamical evolution of the asteroid belt has already 
been the object of a review chapter by Petit et al. (2002) in 
the Asteroid III book. This review has therefore an impor¬ 
tant overlap with that chapter. Nevertheless, both our ob¬ 
servational knowledge of the asteroid belt and our theoreti¬ 
cal understanding of Solar System evolution have improved 
significantly since the early 2000, providing an emerging 
view of a very dynamic early Solar System, in which vari- 
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Fig. 1.— The points show mean proper eccentricity (circles) and mean proper inclination (squares) for the D > 100 km asteroids, 
divided into three bins of semi major axis. The error bars show the l-er standard deviation. There is little systematic difference in 
excitation across the main belt. The slightly increase of inclination from the inner to the outer belt is due to the effect of the g = g 6 
secular resonance (see sect. 3), which most strongly affects high inclination asteroids in the inner belt. 


ous episodes of planet migration played a fundamental role 
in sculpting the small body reservoirs and displacing plan- 
etesimals far from their original birth places. Thus this 
chapter will present in greater details models proposed af¬ 
ter 2002, focusing on their implications for asteroid science. 
Moreover, when reviewing models already presented in Pe¬ 
tit et al., we will refer to numerical simulations of these 
models made after the Petit et al. chapter. 

2. Observational constraints on the primordial evolu¬ 
tion of the asteroid belt 

The observational constraints most useful for recon¬ 
structing the formation and evolution of the asteroid belt are 
those related to large asteroids (larger than ~ 50-100 km in 
diameter). In fact, it has been argued that these asteroids 
are the most likely to be “pristine” in the sense that they 
were not generated in large numbers in collisional break¬ 
up events of larger parent bodies (Bottke et al. 2005a; see 
chapter by Bottke et al.), nor have they been affected by gas 
drag and other non-gravitational forces (e.g., the Yarkovsky 
effect; see chapter by Vokrouhlicky et al.). Moreover, there 
is an emerging view that the first planetesimals were big, 
with a preferred diameter in the range mentioned above 
(Morbidelli et al., 2009; see chapter by Johansen et al.). 
Thus, throughout this chapter we will limit our discussion 
to the properties of large asteroids and refer to smaller as¬ 
teroids only when explicitly mentioned. 

A key major characteristic of the asteroid belt popula¬ 
tion is the orbital excitation, i.e. the fact that the eccen¬ 
tricities and inclinations of many asteroidal orbits are quite 
large (e.g. Petit et al. 2002). The median proper inclina¬ 
tion of D > 100 km asteroids is 11 degrees and the median 


proper eccentricity is 0.145. More importantly, the values 
of eccentricities and inclinations of the largest asteroids are 
considerably dispersed, with the former ranging between 0 
and 0.30, while the latter ranges between 0 and 33 degrees 
(see Fig. 1). It has been shown that asteroids of modest 
inclinations (i < 20°) fill the entire orbital space available 
for long-term dynamical stability, though some stable re¬ 
gions are more densely populated than others (Minton and 
Malhotra, 2009, 2011). The reader should be aware that, 
whatever the preferred formation mechanism (see Johansen 
et al. chapter), planetesimals are expected to have formed 
on circular and co-planar orbits. Thus, one or more dynam¬ 
ical excitation mechanism(s) within the primordial asteroid 
belt were needed to stir up eccentricities and inclinations to 
randomly dispersed values. Asteroid eccentricities and in¬ 
clinations do not show a strong dependence on semi major 
axis (Fig. 1) 

A second fundamental characteristic of the asteroid belt 
is the partial mixing of taxonomic classes. Asteroids can 
be grouped into many taxonomic classes on the basis of 
their visual and infra-red spectroscopic signatures (Tholen, 
1984; Bus and Binzel, 2002; DeMeo et al., 2009). As 
shown first by Gradie and Tedesco (1982) for the largest 
asteroids, the inner belt is dominated by S-complex aster¬ 
oids, many of which are probably related to the meteorites 
known as ordinary chondrites (Binzel et al. 1996; see also 
the chapter by Vernazza et al.). The central belt (2.5-3.2 
AU) is dominated by C-complex asteroids, probably related 
to carbonaceous chondrites (Burbine et al., 2002; see also 
the chapters by DeMeo et al. and Rivkin et al.). The Cy- 
beles asteroids (3.2-3.7 AU), the Hilda asteroids (in the 3/2 
mean motion resonance with Jupiter) and the Jupiter Trojan 
asteroids (in the 1/1 resonance with Jupiter) are dominated 
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Fig. 2.— The relative distribution of large asteroids (D > 50 km) of different taxonomic types as originally observed by Gradie and 
Tedesco (1982). Further works by Mothe-Diniz et al. (2003), Carvano et al. (2010) and DeMeo and Carry (2014) demonstrate that the 
level of mixing increases for smaller asteroid sizes. 


by P-and D-type asteroids (see chapter by Emery et al.). 
The C2 ungrouped meteorite “Tagish Lake” has been pro¬ 
posed to be a fragment of a D-type asteroid (Hiroi et al., 
2001 ). 

This stratification of the main belt makes intuitive sense 
in terms of a general view that proto-planetary disks should 
have temperatures decreasing with increasing distance from 
the central star. In fact, ordinary chondrites are less abun¬ 
dant in organics and water than carbonaceous chondrites 
and therefore are more likely to have formed in a warmer 
part of the disk. The small water content in ordinary chon¬ 
drites, well below the solar proportion, suggests that these 
bodies accreted closer to the Sun than the snowline. The 
fact that some water is nevertheless present is not in con¬ 
tradiction with this statement. A small amount of water 
could have been accreted by collisions with primitive bod¬ 
ies scattered or drifting into the inner part of the disk. At 
the opposite extreme, the Cl meteorites show no chemical 
fractionation relative to the solar composition, except H, C, 
N, O and all noble gases, suggesting that they formed in a 
region of the disk where the temperature was low enough to 
allow the condensation of most elements. 

As shown in Fig. 2, however, asteroids of different taxo¬ 
nomic types are partially mixed in orbital semi major axis, 
which smears the trend relating physical properties to he¬ 
liocentric distance. The mixing of taxonomic type should 
not be interpreted as the existence of asteroids of interme¬ 


diate physical properties between those of adjacent types; 
instead it is due to the coexistence of asteroids of different 
types with various relative proportions at each value of semi 
major axis. Some mixing could come from the fact that 
the thermal and chemical compositional properties of the 
disk probably changed over time. However, given that no 
systematic differences in accretion ages is observed among 
the main group of chondrites (Villeneuve et al., 2009), it 
is more likely that some mechanism, possibly the same 
that excited the orbital eccentricities and inclinations, also 
changed somewhat in a random fashion the original semi 
major axes of the bodies, causing the observed partial mix¬ 
ing. 

The asteroid belt contains overall very little mass. From 
the direct determination of the masses of the largest as¬ 
teroids and an estimate of the total mass of the ring of 
bodies which cannot be individually ’’weighted”, based on 
the collective gravitational perturbations exerted on Mars, 
Krasinsky et al. (2002), Kuchynka and Folkner (2013) and 
Somenzi et al. (2010) concluded that the total mass con¬ 
tained in the asteroid belt is ~ 4.5 x 10 -4 Earth masses 
(Af®). This value is very low compared to that estimated to 
have originally existed in the primordial asteroid belt region 
by interpolating the mass densities required to form the ter¬ 
restrial planets and the core of Jupiter at both ends of the 
belt (Weidenschilling, 1977), which is of the order of 1 
(within a factor of a few). Thus, the mass in the asteroid belt 
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has potentially been depleted by three orders of magnitude 
compared to these expectations. 

We can glean insights into how the primordial belt lost 
its mass by investigating what we know about its collisional 
evolution. The collisional history of asteroids is the subject 
of the chapter by Bottke et al., but we report the highlights 
here that are needed for this discussion. In brief, using 
a number of constraints, Bottke et al. (2005a) concluded 
that the integrated collisional activity of the asteroid belt is 
equivalent to the one that would be produced at the current 
collisional rate over 8-10 Gy. 

This result has several implications. First, it strongly 
suggests that the three orders of magnitude mass deple¬ 
tion could not come purely from collisional erosion; such 
intense comminution would violate numerous constraints. 
Second, it argues that the mass depletion of the asteroid belt 
occurred very early. This is because, once the eccentrici¬ 
ties and inclinations are excited to values comparable to the 
current ones, for a given body every million year spent in 
an asteroid belt 1,000 times more populated brings a num¬ 
ber of collisions equivalent to that suffered in 1 Gy within 
the current population. For this reason, the third implica¬ 
tion is that the dynamical excitation and the mass depletion 
event almost certainly coincided. This argues that the real 
dynamical excitation event was stronger than suggested by 
the current distribution of asteroid eccentricities. One way 
to reconcile a massive asteroid belt with this scenario is to 
assume that over 99% of the asteroids had their orbits so ex¬ 
cited that they left the asteroid belt forever (hence the mass 
depletion). This would make the eccentricities (and to a 
lesser extent the inclinations) we see today to be those de¬ 
fined by the the lucky survivors, namely the bodies whose 
orbits were excited the least. 

Using these constraints, we discuss in the next session 
the various models that have been proposed for the primor¬ 
dial sculpting of the asteroid belt. 

3. Models of primordial evolution of the asteroid belt 
3.1. Early models 

The first attempts to explain the primordial dynamical 
excitation of the asteroid belt were made by Heppenheimer 
(1980) and Ward (1981) who proposed that secular reso¬ 
nances swept through the asteroid belt region during the 
dissipation of gas in the proto-planetary disk. Secular res¬ 
onances occur when the precession rate of the orbit of an 
asteroid is equal to one of the fundamental frequencies of 
precession of the orbits of the planets. There are two angles 
that characterize the orientation of an orbit in space, the lon¬ 
gitude of perihelion ( vj ) and the longitude of the ascending 
node (12), each of which can precess at different rates de¬ 
pending on the gravitational effects of the other planets and 
nebular gas (if present). The resonances that occur when 
the precession rates of the longitudes of perihelion of an 
asteroid (denoted by g) and of a planet are equal to each 
other excite the asteroid’s eccentricity. Similarly, the res¬ 
onances occurring when the precession rates of the longi¬ 


tudes of node of an asteroid (denoted by s) and of a planet 
are equal to each other excite the asteroid’s inclination. In 
the case of asteroids in the main belt, the planets’ precession 
frequencies that most influence their dynamics are those as¬ 
sociated with the orbits of lupiter and Saturn. These are 
called and g§ for the longitude of perihelion precession 
(the former dominating in the precession of the perihelion 
of lupiter, the latter in that of Saturn), and s e for the longi¬ 
tude of the node precession (both the nodes of Jupiter and 
Saturn precess at the same rate, if defined relative to the in¬ 
variable plane, defined as the plane orthogonal to their total 
angular momentum vector). 

The dissipation of gas from the proto-planetary disk 
changes the gravitational potentials that the asteroids and 
planets feel, and hence changes the precession rates of their 
orbits. Given that the planets and asteroids are at differ¬ 
ent locations, they will be affected somewhat differently by 
this change of gravitational potential and consequently their 
precession rates will not change proportionally. It is there¬ 
fore possible that secular resonances sweep through the as¬ 
teroid belt as the gas dissipates. This means that every as¬ 
teroid, whatever its location in the belt, first has orbital pre¬ 
cession rates slower than the g -,, g& frequencies of Jupiter 
and Saturn when there is a lot of gas in the disk, then en¬ 
ters resonance (g = g 5 or g = g lt ) when some appropriate 
fraction of the gas has been removed, and eventually is no 
longer in resonance (its orbital precession frequency being 
faster than those of the giant planets: i.e. g > go) after all 
the gas has disappeared. The same occurs for the asteroid’s 
nodal frequency s relative to the planetary frequency Sq. 
This sweeping of perihelion and nodal secular resonances 
has the potential to excite the orbital eccentricities and in¬ 
clinations of all asteroids. 

This mechanism of asteroid excitation due to disk dis¬ 
sipation has been revisited with numerical simulations in 
Lemaitre and Dubru (1991), Lecar and Franklin (1997), Na- 
gasawa et al. (2000,2001,2002), Petit et al. (2002), and 
finally by O’Brien et al. (2007). Nagasawa et al. (2000) 
found that of all the scenarios for gas depletion they stud¬ 
ied (uniform depletion, inside-out, and outside-in), inside- 
out depletion of the nebula was most effective at exciting 
eccentricities and inclinations of asteroids throughout the 
main belt. However, they (unrealistically) assumed that the 
nebula coincided with the ecliptic plane. Proto-planetary 
disks can be warped, but they are typically aligned with the 
orbit of the locally dominant planet (Mouillet et al., 1997). 
Thus, there is no reason that the gaseous disk in the aster¬ 
oid belt region was aligned with the current orbital planet 
of the Earth (which was not formed yet). Almost certainly 
it was aligned with the orbits of the giant planets. Taking 
the invariable plane (the plane orthogonal to the total an¬ 
gular momentum of the Solar System) as a proxy of the 
original orbital plane of Jupiter and Saturn, Nagasawa et 
al. (2001, 2002) found that the excitation of inclinations 
would be greatly diminished. Furthermore, since nebular 
gas in the inside-out depletion scenario would be removed 
from the asteroid belt region before the resonances swept 
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through it, there would be no gas drag effect to help deplete 
material from the main belt region. 

The work of O’Brien et al. (2007) accounted for the fact 
that the giant planets should have had orbits significantly 
less inclined and eccentric than their current values when 
they were still embedded in the disk of gas, because of the 
strong damping that gas exerts on planets (Cresswell et al. 
2008; Kley and Nelson, 2012). They concluded that sec¬ 
ular resonance sweeping is effective at exciting eccentrici¬ 
ties and inclinations to their current values only if gas is re¬ 
moved from the inside-out and very slowly, on a timescale 
of ~20 My. This gas-removal mode is very different from 
our current understanding of the photo-evaporation process 
(Alexander et al., 2014), and inconsistent with observations 
suggesting that disks around solar-type stars have lifetimes 
of only 1-10 My, with an average of ~3 My (eg. Strom et 
al. 1993; Zuckerman et al. 1995; Kenyon and Hartmann 
1995; Haisch et al. 2001). 

Earlier studies found that the final eccentricities of the 
asteroids are quite randomized because two perihelion sec¬ 
ular resonances sweep the entire asteroid belt in sequence- 
first the resonance g = g-„ then the resonance g = g§. 
The first resonance excites the eccentricities of the asteroids 
from zero to approximately the same value, but the second 
resonance, sweeping an already excited belt, can increase 
or decrease the eccentricity depending on the position of 
the perihelion of each asteroid at the time of the encounter 
with the resonance (Ward et al., 1976; Minton and Malho- 
tra, 2011). O’Brien et al. (2007) found that when Jupiter 
and Saturn were on orbits initially closer together, as pre¬ 
dicted by the Nice Model (e.g., Tsiganis et al. 2005), the 
resonance with frequency g<, would only sweep part of the 
outer belt, leading to less randomization of eccentricities 
in the inner belt. In all studies in which the mid-plane of 
the proto-planetary disk of gas coincides with the invariable 
plane of the solar system find that the final inclinations tend 
to have comparable values. This is because there is only one 
dominant frequency (se) in the precession of the nodes of 
Jupiter and Saturn and hence there is only one nodal secu¬ 
lar resonance and no randomization of the final inclinations 
of the asteroids. Clearly, this is in contrast with the obser¬ 
vations. For all these problems, the model of secular reso¬ 
nance sweeping during gas removal is no longer considered 
to be able to explain, alone, the excitation and depletion of 
the primordial asteroid belt. 

An alternative model for the dynamical excitation of the 
asteroid belt was proposed by Ip (1987). In this model, 
putative planetary embryos are scattered out of the Jupiter 
region and cross the asteroid belt for some timescale be¬ 
fore being ultimately dynamically ejected from the So¬ 
lar System. If the embryos are massive enough, their re¬ 
peated crossing of the asteroid belt can excite and ran¬ 
domize the eccentricities and inclinations of the asteroids, 
through close encounters and secular effects. That scenario 
has been revisited by Petit et al. (1999), who found that, 
whatever the mass of the putative embryos, the resulting 
excitation in the asteroid belt ought to be very unbalanced. 


Excitation would be much stronger in the outer belt than 
in the inner belt (because the embryos come from Jupiter’s 
region) and it would be much stronger in eccentricity than 
in inclination. By contrast, the main asteroid belt shows no 
such trend (see Fig. 1). So, again, this model has since 
been abandoned. If massive embryos have been scattered 
from Jupiter’s zone, they must have crossed the asteroid belt 
very briefly so that their limited effects could be completely 
overprinted by other processes, such as those discussed be¬ 
low. 

3.2. Wetherill’s model 

The first comprehensive model of asteroid belt sculpt¬ 
ing, which linked the evolution of the asteroid belt with the 
process of terrestrial planet formation, was that proposed 
by Wetherill (1992) and later simulated in a number of sub¬ 
sequent papers (e.g.. Chambers and Wetherill, 1998; Pe¬ 
tit et al., 2001, 2002; O’Brien et al., 2006, 2007). In this 
model, at the time gas was removed from the system, the 
proto-planetary disk interior to Jupiter consisted of a bi- 
modal population of planetesimals and planetary embryos, 
the latter with masses comparable to those of the Moon 
or Mars. Numerical simulations show that, under the ef¬ 
fect of the mutual perturbations among the embryos and the 
resonant perturbations from Jupiter, embryos are generally 
cleared from the asteroid belt region, whereas embryos col¬ 
lide with each other and build terrestrial planets inside of 
2 AU. While they are still crossing the asteroid belt, the 
embryos also excite and eject most of the original resident 
planetesimals. Only a minority of the planetesimals (and 
often no embryos) remain in the belt at the end of the ter¬ 
restrial planets formation process, which explains the mass 
depletion of the current asteroid population. The eccentric¬ 
ities and inclinations of the surviving asteroids are excited 
and randomized, and the remaining asteroids have gener¬ 
ally been scattered somewhat relative to their original semi 
major axes. A series of simulation snapshots demonstrating 
this process is shown in Figure 3. 

Whereas earlier simulations assumed that Jupiter and 
Saturn were originally on their current orbits, O’Brien et al. 
(2006, 2007) performed simulations with Jupiter and Sat¬ 
urn on the low-inclination, nearly circular orbits predicted 
in the Nice Model. The resulting asteroids from a set of 
simulations with these initial conditions are shown in Fig¬ 
ure 4. Overall, the range of values compare well with those 
observed for the real asteroids, although the final inclination 
distribution is skewed towards large inclinations. The rea¬ 
son for this is that it is easier to excite a low-inclination as¬ 
teroid to large eccentricity and remove it from the belt than 
it is for a high-inclination asteroid, because the encounter 
velocities with the embryos are slower and more effective 
in deflecting the low-inclination asteroid’s orbit. Also, with 
the giant planets on nearly circular orbits, it takes longer to 
clear embryos from the asteroid belt, allowing more time to 
excite asteroids to large inclinations. 

As noted earlier, the surviving asteroids have their or- 
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Fig. 3.— Snapshots of the evolution of the solar system and of the asteroid belt in a simulation of WetherilFs model performed in 
O'Brien et al. (2006) and assuming Jupiter and Saturn on initial quasi-circular orbits. Each panel depicts the eccentricity vs. semi 
major axis distribution of the particles in the system at different times, labeled on top. Planetesimals are represented with gray dots and 
planetary embryos by black circles, whose size is proportional to the cubic root of their mass. The solid lines show the approximate 
boundaries of the current main belt. 


bital semi major axes displaced from their original values, 
as a result of the embryos’ gravitational scattering. O’Brien 
et al. (2007) found that the typical change in semi major 
axis is of the order of 0.5 AU (comparable to earlier simula¬ 
tions), with a tail extending to 1 - 2 AU, and the semi major 
axis can be either decreased or increased. This process can 
explain the partial mixing of taxonomic types. As shown in 
Fig. 2 the distribution of the S-type and C-type asteroids has 
a Gaussian-like shape, with a characteristic width of ~0.5 
AU. Thus, if one postulates that all S-type originated from 
the vicinity of 2 AU and all C-type originated in the vicinity 
of 3 AU, Wetherill’s model explains the current distribution. 

3.3. The Grand Tack model 

A more recent, alternative model to WetherilFs is the 
so-called Grand Tack scenario, proposed in Walsh et al. 
(2011). Initially the Grand Tack scenario had not been de¬ 
veloped to explain the asteroid belt, but to answer two ques¬ 


tions left open by WetherilFs model: why is Mars so small 
relative to the Earth? Why is Jupiter so far from the Sun de¬ 
spite planets having a tendency to migrate inwards in proto¬ 
planetary disks? Nevertheless, this scenario has profound 
implications for the asteroid belt, as we discuss below. 

The Grand Tack scenario is built on results from hy¬ 
drodynamics simulations finding that Jupiter migrates to¬ 
wards the Sun if it is alone in the gas-disk, while it mi¬ 
grates outward if paired with Saturn (Masset and Snell- 
grove, 2001; Morbidelli and Crida, 2007; Pierens and Nel¬ 
son, 2008; Pierens and Raymond, 2011; D’angelo and 
Marzari, 2012). Thus, the scenario postulates that Jupiter 
formed first. As long as the planet was basically alone, Sat¬ 
urn being too small to influence its dynamics, Jupiter mi¬ 
grated inwards from its initial position (poorly constrained 
but estimated at ~3.5 AU) down to 1.5 AU. Then, when Sat¬ 
urn reached a mass close to its current one and an orbit close 
to that of Jupiter, Jupiter reversed migration direction (aka 
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Fig. 4.— The final eccentricities and inclinations of asteroids in Wetherill’s (1992) model (black dots), according to the simulations 
presented in O'Brien et al. (2007). For comparison, the observed distribution of large asteroids is depicted with gray dots. 


it ’’tacked”, hence the name of the model) and the pair of 
planets started to move outwards. This outward migration 
continued until the final removal of gas in the disk, which 
the model assumes happened when Jupiter reached a dis¬ 
tance of ~5.5 AU. The migration of the cores of giant plan¬ 
ets is still not fully understood (see Kley and Nelson, 2012 
for a review). Thus, the Grand Tack model comes in two 
flavors. In one, Saturn, while growing, migrates inwards 
with Jupiter. In another, Saturn is stranded at a no-migration 
orbital radius until its mass exceeds 50 Mj, (Bitsch et al., 
2014); then it starts migrating inwards and it catches Jupiter 
in resonance because it migrates faster. Both versions exist 
with and without Uranus and Neptune. All these variants 
are described in Walsh et al. (2011); the results are very 
similar in all these cases, which shows the robustness of 
the model, at least withing the range of tested possibilities. 
The scheme presented in Fig. 5. has been developed in the 
framework of the first “flavor”. 

Assuming that Jupiter formed at the snowline (a usual 
assumption to justify the large mass of its core and its fast 
formation), the planetesimals that formed inside its initial 
orbit should have been mostly dry. It is therefore reason¬ 
able to associate these planetesimals (whose distribution is 
sketched as a dashed area in Fig. 5) with the S-type as¬ 
teroids and other even dryer bodies (enstatite-type. Earth 
precursors etc.). During its inward migration, Jupiter pene¬ 
trates into the disk of these planetesimals. In doing so, most 
planetesimals (and planetary embryos) are captured in mean 
motion resonances with Jupiter and are pushed inwards, in¬ 
creasing the mass density of the inner part of the disk. How¬ 
ever, some 10% of the planetesimals are kicked outwards by 
an encounter with Jupiter, reaching orbits located beyond 
Saturn, which collectively have an orbital (a,e) distribution 
that is typical of a scattered disk (i.e. with mean eccentricity 
increasing with semi major axis). In semi major axis range, 
this scattered disk overlaps with the inner part of the disk of 


primitive bodies (whose distribution is sketched as a dotted 
area in Fig. 5), which are initially on circular orbits beyond 
the orbit of Saturn. These bodies, being formed beyond the 
snowline, should be rich in water ice and other volatile ele¬ 
ments, and therefore it is again reasonable to associate them 
with C-type asteroids. 

After reaching ~ 1.5 AU (this value is constrained by the 
requirement to form a small Mars and a big Earth; Walsh et 
al., 2011; Jacobson et al., 2014; Jacobson and Morbidelli, 
2014), Jupiter reverses its migration direction and begins 
its outward migration phase, during which the giant planets 
encounter the scattered S-type disk, and then also the prim¬ 
itive C-type disk. Some of the bodies in both populations 
are scattered inwards, reach the asteroid belt region and are 
implanted there as Jupiter moves out of it. 

The final orbits of the planetesimals, at the end of the 
outward migration phase, are shown in Fig. 6. A larger 
dot size is used to highlight the planetesimals trapped in 
the asteroid belt region and distinguish them from those in 
the inner solar system or at too large eccentricity to be in 
the asteroid belt. Notice that the asteroid belt is not empty, 
although it has been strongly depleted (by a factor of sev¬ 
eral hundred relative to its initial population). This result 
is not trivial. One could have expected that Jupiter migrat¬ 
ing through the asteroid belt twice (first inwards then out¬ 
wards) would have completely removed the asteroid popu¬ 
lation, invalidating the Grand Tack scenario. The eccentric¬ 
ities and the inclinations of the particles in the asteroid belt 
are excited and randomized. The S-type particles (black) 
are found predominantly in the inner part of the belt and the 
C-type particles (gray) in the outer part, but there is a wide 
overlapping region where both are present. This is qualita¬ 
tively consistent with what is observed. 

As discussed above, the Grand Tack scenario solves 
open problems in Wetherill’s model. The small mass of 
Mars is explained as a result of the disk of the remaining 
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Fig. 5.— A scheme showing the Grand Tack evolution of Jupiter and Saturn and its effects on the asteroid belt. The three panels show 
three evolutionary states, in temporal sequence. First the planet migrate inwards then, when Saturn reaches its current mass, they move 
outwards. The dashed and dotted areas schematize the (a,e) distributions of S-type and C-type asteroids respectively. The dashed and 
dotted arrows in the lower panel illustrate the injection of scattered S-type and C-type asteroids into the asteroid belt during the final 
phase of outward migration of the planets. 


solid material being truncated at ~1 AU (Hansen, 2009; 
Walsh et al., 2011). Fischer and Ciesla (2014) reported 
that they could obtain a small-mass Mars in a few percent 
of simulations conducted in the framework of Wetherill’s 
model. However, the rest of the planetary system in these 
simulations does not resemble the real terrestrial planet sys¬ 
tem (Jacobson and Walsh, 2015). For instance another mas¬ 
sive planet is formed in Mars-region or beyond. The out¬ 
ward migration of Jupiter explains why the giant planets in 
our solar system are so far from the Sun, whereas most gi¬ 
ant planets found so far around other stars are located at 1-2 
AU. For all these reasons, one can consider the Grand Tack 
model more as an improvement of Wetherill’s model than 
an alternative, because it is built in the same spirit of linking 
the asteroid belt sculpting to the evolution of the rest of the 
solar system (terrestrial planet formation, giant planet mi¬ 
gration - the latter being still unknown at Wetherill’s time). 

It is nevertheless interesting to compare the Grand Tack 
model and Wetherill’s model on the basis of the final as¬ 
teroid belts that they produce. Comparing Fig. 6 with Fig. 
4, it is apparent that the Grand Tack model provides a bet¬ 
ter inclination distribution, more uniform that Wetherill’s, 
but it produces a worse eccentricity distribution, which is 
now more skewed towards the upper eccentricity boundary 
of the asteroid belt. 

As we will see in Sect. 3.5 however, the eccentricity dis¬ 
tribution can be remodeled somewhat during a later evolu¬ 


tionary phase of the solar system. This is also partially true 
also for the inclination distribution. So, for what concerns 
the eccentricity and inclination distributions one might de¬ 
clare a tie in the competition between the two models. 

The Grand Tack model makes it conceptually easier to 
understand the significant differences between S-type and 
C-type asteroids and their respective presumed daughter 
populations: the ordinary and carbonaceous chondrites. In 
fact, in the Grand Tack model these two populations are 
sourced from clearly distinct reservoirs on either sides of 
the snowline. Instead, in Wetherill’s model these bodies 
would have formed just at the two ends of the asteroid belt, 
so less than 1 AU apart. Despite such a vast difference in 
predicted formation locations for these two populations the 
debate is open. Some authors (e.g. Alexander et al., 2012) 
think that bodies formed in the giant planet region would 
be much more similar to comets than to asteroids, others 
(Gounelle et al., 2008) argue that there is a continuum be¬ 
tween C-type asteroids and comets and a clear cleavage 
of physical properties between ordinary and carbonaceous 
chondrites. We review the available cosmochemical con¬ 
straints and their uncertain compatibility with the model in 
sect. 13.41 

A clear distinction between the Grand Tack model and 
Wetherill’s model is that the former provides a faster and 
more drastic depletion of the asteroid belt. This point is 
illustrated in Fig. 7, showing the fraction of the initial as- 
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Fig. 6.— Final semi major axis, eccentricity and inclination distribution of bodies surviving the inward and outward migration of 
Jupiter and Saturn. The black particles were originally placed inside of the initial orbit of Jupiter and the gray particles outside of the 
initial orbit of Saturn. The particles finally trapped in the asteroid belt are depicted with larger symbols than the others. The dashed curve 
in the lower panel shows the approximate boundaries of the asteroid belt inward of the 2/1 resonance with Jupiter. This final distribution 
was achieved in the simulations of Walsh et al. (2011) accounting only for Jupiter and Saturn (i.e., not including Uranus and Neptune) 
moving together in the 2/3 resonance, as sketched in Fig. 5. 


teroid population that is in the main belt region at any time. 

The Grand Tack scenario depletes the asteroid belt down 
to 0.3%, and does so basically in 0.1 My. Assuming that the 
final asteroid belt consisted of one current asteroid belt mass 
in S-type asteroids and three current asteroid belt masses 
in C-type asteroids (the reason for 4x more total mass in 
the asteroid belt will be clarified in Sect. |3.6| ), this implies 
that the asteroid belt at t=0 should have contained 0.6 Earth 
masses in planetesimals (the rest in embryos). Also, a cal¬ 
culation of the collision probability of the asteroids as a 
function of time (both among each other and with the plan¬ 
etesimals outside of the asteroid belt) shows that the inte¬ 
grated collisional activity suffered by the surviving aster¬ 
oids during the first 200 My would not exceed the equiva¬ 
lent of 4 Gy in the current population. Thus, assuming that 
the exceeding factor of 4 in the asteroid population is lost 
within the next 500 My (see Sects. |3.5| and |3.6| ), the inte¬ 
grated collisional activity of asteroids throughout the entire 
solar system age would probably remain within the 10 Gy 
constraint described in Section 2. In contrast, Wetherill’s 


model depletes the asteroid belt on a timescale of 100 My. 
Also, about 2-3 % of the initial population remains in the 
belt at the end. Thus, to be consistent with constraints on 
the current population and its integrated collisional activity, 
the initial mass in planetesimals in the asteroid belt region 
should have been no larger than 200 times the current as¬ 
teroid belt mass, or less than one Mars mass (Bottke et al., 
2005b). 

3.4. Are cosmochemical constraints consistent with the 
Grand Tack model? 

The Grand Tack model predicts that C-type asteroids 
have been implanted into the asteroid belt from the giant 
planets region. Is this supported or refuted by cosmochem¬ 
ical evidence? 

Although there is a spread in values, the D/H ratios of 
carbonaceous chondrites (with the exception of CR chon¬ 
drites) are a good match to Earth’s water (Alexander et al., 
2012). Oort cloud comets are usually considered to have 
formed in the giant planet region (e.g. Dones et al., 2004). 
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Fig. 7.— The depletion of the asteroid belt in Wetherill’s model (black) and Grand Tack model (gray). In the Grand Tack model notice 
the bump between 0.1 and 0.6 My due to the implantation of primitive objects into the main belt. Overall, the depletion of the asteroid 
belt is faster and stronger in the Grand Tack model. 


The D/H ratio was measured for the water from seven Oort 
cloud comets see Bockelee-Morvan et al., 2012 and refer¬ 
ences in that paper). All but one (comet 153P/Ikeya-Zhang; 
Biver et al., 2006) have water D/H ratios of about twice 
higher than chondritic. This prompted Yang et al. (2013) 
to develop a model where the D/H ratio of ice in the gi¬ 
ant planet region is high. However, Brasser et al. (2007) 
showed that comet-size bodies could not be scattered from 
the giant planet region into the Oort cloud in the presence of 
gas drag (i.e., when the giant planets formed), and Brasser 
and Morbidelli (2013) demonstrated that the Oort cloud 
population is consistent with an origin from the primordial 
trans-Neptunian disk at a later time. The recent measure¬ 
ment (Altwegg et al., 2014) of a high D/H ratio for the ice 
of comet 67P/Tchourioumov-Guerassimenko, which comes 
from the Kuiper belt, supports this conclusion by showing 
that there is no systematic difference between Oort cloud 
comets and Kuiper belt comets. So, care should be taken 
in using Oort cloud comets as indicators of the D/H ratio 
in the giant planet region. Conflicting indications on the 
local D/H ratio come from the analysis of Saturn’s moons. 
Enceladus’ D/H ratio is roughly twice Earth’s (Waite et al. 
2009), but Titan’s D/H ratio is Earth-like (Coustenis et al. 


2008; Abbas et al. 2010; Nixon et al. 2012). 

Alexander et al. also noticed a correlation between D/H 
and C/H in meteorites and interpreted it as evidence for an 
isotopic exchange between pristine ice and organic matter 
within the parent bodies of carbonaceous chondrites. From 
this consideration, they argued that the original water reser¬ 
voir of carbonaceous asteroids had a D/H ratio lower than 
Titan, Enceladus or any comet, again making asteroids dis¬ 
tinct from bodies formed in the giant planet region and be¬ 
yond. However, a reservoir of pristine ice has never been 
observed; the fact that Earth’s water and other volatiles are 
in chondritic proportion (Marty, 2012; however see Halli- 
day, 2013) means that carbonaceous chondrites -wherever 
they formed- reached their current D/H ratios very quickly, 
before delivering volatiles to Earth. Possibly, also the D/H 
ratio measured for comets and satellites might have been 
the result of a similar rapid exchange between a pristine ice 
and the organic matter. 

Another isotopic constraint comes from the Nitrogen 
isotope ratio. Comets seem to have a rather uniform 
15 N/ 14 N (Rousselot et al., 2014). Even the comets with 
a chondritic D/H ratio (e.g. Hartley 2; Hartogh et al. 2011) 
have a non-chondritic 15 N/ 14 N ratio (Meech et al., 2011). 
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The 15 N/ 14 N ratio, however, is only measured in HCN or 
NH 2 , never in molecular Nitrogen N 2 . Titan has a cometary 
15 N/ 14 N as well (in this case measured in N 2 ; Niemann et 
al., 2010). Here, again, a few caveats are in order. First, it 
is difficult to relate the composition of a satellite, born from 
a circum-planetary disk with its own thermal and chemical 
evolution, to the composition of bodies born at the same 
solar distance but on heliocentric orbits. Second, it is un¬ 
clear whether any comets for which isotope ratios have been 
measured originate from the giant planet region, as opposed 
to the trans-Neptunian disk (Brasser and Morbidelli, 2013). 
We remind that hot-spots with large 15 N/ 14 N ratios are also 
found in primitive meteorites (Busemann et al., 2006). 

Arguments in favor of an isotopic similarity between car¬ 
bonaceous chondrites and comets come from the analysis 
of micro-meteorties. Most of micro-meteorites (particles of 
about ~ 100/rm collected in the Antarctic ice) have chon- 
dritic isotopic ratios for H, C, N, O (with the exceptions 
of the ultra-carbonaceous ones, which have a large D/H ra¬ 
tio, comparable to that in the organics of some chondrites, 
but which constitute only a very minority of micromete¬ 
orite population - Duprat et al., 2010). Yet, according to 
the best available model of the dust distribution in the inner 
solar system (Nesvorny et al., 2010), which is compelling 
given that it fits the zodiacal light observations almost per¬ 
fectly, most of the dust accreted by the Earth should be 
cometary, even when the entry velocity bias is taken into 
account. Similarly, from orbital considerations, Gounelle 
et al. (2006) concluded that the Cl meteorite Orgueil is 
a piece of a comet. Also compelling for a continuum be¬ 
tween chondrites and comets comes from the examination 
of comet Wild 2 particles brought to Earth by the Stardust 
Mission - e.g., Zolensky et al. 2006, 2008.; Ishii et al. 2008; 
Nakamura et al., 2008 

These considerations suggest that if one looks at their 
rocky components, comets and carbonaceous asteroids are 
very similar from the composition and isotopic point of 
view, if not indistinguishable. 

Finally, it has been argued that, if the parent bodies 
of carbonaceous chondrites had accreted among the giant 
planets they would have contained ~ 50% water by mass. 
Instead, the limited hydrous alteration in carbonaceous me¬ 
teorites suggests that only about 10% of the original mass 
was in water (Krot, 2014; however see Alexander et al., 
2010). However, a body’s original water content cannot 
easily be estimated from its aqueous alteration. Even if al¬ 
teration is complete, there is a finite amount of water that 
the clays can hold in their structures. Thus, the carbona¬ 
ceous parent bodies may have been more water-rich than 
their alteration seems to imply. In fact, the discoveries of 
main belt comets releasing dust at each perihelion passage 
(Hsieh and Jewitt, 2006), of water ice on asteroids Themis 
(Campins et al., 2010; Rivkin and Emery, 2010) and Cy- 
bele (Hargrove et al., 2012), and of vapor plumes on Ceres 
(Kuppers et al., 2014), show that C-type asteroids are more 
rich in water than their meteorite counterpart seems to sug¬ 
gest, supporting the idea that they might have formed near 


or beyond the snowline. Meteorites may simply represent 
rocky fragments of bodies that were far wetter/icier. 

Clearly, the debate on whether carbonaceous asteroids 
really come from the giant planet region as predicted by 
the Grand Tack model is wide open. More data are needed 
from a broader population of comets. The investigation of 
main belt comets, both remote and in-situ, and the Dawn 
mission at Ceres will be key to elucidate the real ice-content 
of carbonaceous asteroids and their relationship with classic 
’’comets”. 

3.5. The Nice model: a second phase of excitation and 
depletion for the asteroid belt 

Fig. 7. seems to suggest that after ~100 My the aster¬ 
oid belt had basically reached its final state. However, at 
this time the orbits of the giant planets were probably still 
not the current ones. In fact, the giant planets are expected 
to have emerged from the gas-disk phase in a compact and 
resonant configuration as a result of their migration in the 
gas-dominated disk. This is true not only in the Grand Tack 
model, but in any model where Jupiter is refrained from mi¬ 
grating inside ^5 AU by whatever mechanism (Morbidelli 
et al., 2007). 

The transition of the giant planets from this early con¬ 
figuration to the current configuration is expected to have 
happened via an orbital instability, driven by the interaction 
of the planets with a massive disk of planetesimals located 
from a few AUs beyond the original orbit of Neptune to 
about 30 AU (Morbidelli et al., 2007; Batygin and Brown, 
2010; Levison et al., 2011; Nesvorny, 2011; Batygin et al., 
2012; Nesvorny and Morbidelli, 2012; see also a precursor 
work by Thommes et al., 1999). In essence, the planetes¬ 
imals disturbed the orbits of the giant planets and, as soon 
as two planets fell off resonance, the entire system became 
unstable. In the simulations, the instability can occur early 
(e.g. Morbidelli et al., 2007) or late (Levison et al., 2011) 
depending on the location of the inner edge of the trans- 
Neptunian disk. 

Constraints suggest that in the real Solar System the in¬ 
stability occurred relatively late, probably around 4.1 Gy 
ago (namely 450 My after gas removal). These constraints 
come primarily from the Moon. Dating lunar impact basins 
is difficult, because it is not clear which samples are re¬ 
lated to which basin (e.g., Norman and Nemchin, 2014). 
Nevertheless, it is clear that several impact basins, proba¬ 
bly a dozen, formed in the 4.1-3.8 Gy period (see Fassett 
and Minton, 2013, for a review). Numerical tests demon¬ 
strate that these late basins (even just Imbrium and Ori- 
entale, whose young ages are undisputed) are unlikely to 
have been produced by a declining population of planetesi¬ 
mals, left-over from the terrestrial planet accretion process, 
because of their short dynamical and collisional lifetimes 
(Bottke et al. 2007). There is also a surge in lunar rock im¬ 
pact ages ~4 Gy ago, which contrasts with a paucity of im¬ 
pact ages between 4.4 and 4.2 Gy (Cohen et al., 2005). This 
is difficult to explain if the bombardment had been caused 
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Fig. 8.— A comparison between the final (a,i) distribution of asteroids if Jupiter and Saturn migrate slowly away from each other (left) 
or jump due to them scattering an ice giant planet (right). In the first case the orbital distribution is inconsistent with that observed, while 
in the second case it is. From Morbidelli et al. (2010). 


by a population of left-over planetesimals slowly declining 
over time. The situation is very similar for the bombard¬ 
ment of asteroids, with meteorites showing a surge in im¬ 
pact ages 4.1 Gy ago and a paucity of ages between 4.2-4.4 
Gy (Marchi et al., 2013). Meteorites also show many im¬ 
pact ages near 4.5 Gy ago, demonstrating that the apparent 
lack of events in the 4.2-4.4 Gy interval is not due to clock 
resetting processes. All these constraints strongly suggest 
the appearance of a new generation of projectiles in the in¬ 
ner solar system about 4.1 Gy ago, which argues that either 
a very big asteroid broke up at that time (Cuk, 2012; but 
such a break-up is very unlikely from collision probabil¬ 
ity arguments and we don’t see any remnant asteroid family 
supporting this hypothesis), or that the dynamical instability 
of the giant planets occurred at that time, partially destabi¬ 
lizing small body reservoirs that had remained stable until 
then. 

Other constraints pointing to the late instability of the gi¬ 
ant planets come from the outer Solar System. If the planets 
had become unstable at the disappearance of the gas in the 
disk, presumably the Sun would still have been in a stellar 
cluster and consequently the Oort cloud would have formed 
more tightly bound to the Sun than it is thought to be from 
the orbital distribution of long period comets (Brasser et al., 
2008, 2012). Also, the impact basins on Iapetus (a satel¬ 
lite of Saturn) have topographies that have relaxed by 25% 
or less, which argues that they formed in a very viscous 
lithosphere; according to models of the thermal evolution of 
the satellite, these basins can not have formed earlier than 
200 My after the beginning of the Solar System (Robuchon 
et al., 2011). 

For all these reasons, it is appropriate to discuss the con¬ 
sequences of the giant planet instability on the asteroid belt, 
after the events described by the Grand Tack or Wetherill’s 
models. In fact, it is important to realize that the model 
of giant planets’ instability (often called the Nice model) is 
not alternative to the models described before on the early 
evolution of the asteroid belt; instead it is a model of the 


subsequent evolution. 

The phase of giant planet instability is very chaotic and 
therefore a variety of orbital evolutions are possible. Nev¬ 
ertheless, the planetary evolutions can be grouped in two 
categories. In the first category, Saturn, Uranus and Nep¬ 
tune have close encounters with each other, but none of 
them have encounters with Jupiter. Saturn typically scat¬ 
ters either Uranus or Neptune outwards and thus it recoils 
inwards. As a result, Uranus and Neptune acquire large ec¬ 
centricity orbits that cross the trans-Neptunian disk. The 
dispersal of the planetesimal disk damps the eccentricities 
of the planets by dynamical friction and drives the plan¬ 
ets’ divergent migration in semi major axis (Tsiganis et al., 
2005). Thus, the planets achieve stable orbits that are well- 
separated from each other. The orbital separation between 
Jupiter and Saturn first decreases, when Saturn recoils, and 
then increases due to planetesimal-driven migration. The 
timescale for the latter is typically ~ lOMy. The slow sep¬ 
aration between the two major planets of the Solar System 
drives a slow migration of secular resonances across the as¬ 
teroid belt (Minton and Malhotra, 2009, 2011) and the ter¬ 
restrial planet region (Brasser et al., 2009). The problem is 
that the resulting orbital distribution in the asteroid belt is 
very uneven, as shown in Fig. 8a, with most asteroids sur¬ 
viving at large inclination (Morbidelli et al., 2010) and the 
orbits of the terrestrial planets become too excited (Brasser 
et al., 2009; Agnor and Lin, 2012). 

In the second category of evolutions, Saturn scatters an 
ice giant planet (Uranus, Neptune or a putative fifth planet) 
inwards, thus recoiling outwards, and then Jupiter scatters 
the ice giant outwards, thus recoiling inwards. The inter¬ 
action with the planetesimals eventually damp and stabilize 
the orbits of the planets. In this evolution, dubbed ’jumping- 
Jupiter’, the orbital separation between Jupiter and Saturn 
initially jumps, when Saturn recoils outwards and Jupiter 
inwards; then there is a final smooth phase of separation, 
due to planetesimal driven migration. In the jump, the secu¬ 
lar resonances can jump across the asteroid belt (Morbidelli 
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et al., 2010) and across the terrestrial planet region (Brasser 
et al., 2009), without disrupting their orbital structure (see 
Fig. 8b). 

The jumping-Jupiter evolution also explains the capture 
of the irregular satellites of Jupiter with an orbital distribu¬ 
tion similar to those of the irregular satellites of the other 
giant planets (Nesvorny et al., 2007, 2014). It can also ex¬ 
plain the capture of Jupiter’s Trojans in uneven proportions 
around the L4 and L5 Lagrangian points (Nesvorny et al., 
2013; see chapter by Emery et al.).So far, no other model 
is capable to achieve these results. For all these reasons, 
nowadays simulated Solar System evolutions are required 
to show a jumping-Jupiter evolution to be declared success¬ 
ful (Nesvorny and Morbidelli, 2012). 

Although it avoids secular resonances sweeping across 
the asteroid belt, a jumping-Jupiter evolution is not without 
consequences for the asteroids. The sudden change in the 
eccentricity of Jupiter (from an initial basically circular or¬ 
bit, like that observed in hydrodynamical simulations of the 
four giant planets evolving in the gaseous proto-planetary 
disk, to one with current eccentricity) changes the forced 
eccentricity felt by the asteroids in their secular evolution. 
Consequently, the proper eccentricities of the asteroids are 
changed. Depending on the value of the longitude of peri¬ 
helion when the forced eccentricity changes, the proper ec¬ 
centricity of an asteroid can decrease or increase. Roughly, 
50% of the asteroids are kicked to larger eccentricities and 
therefore are removed from the asteroid belt. The remain¬ 
ing 50% of the asteroids have their eccentricities reduced. 
This can potentially reconcile the eccentricity distribution 
of asteroids at the end of the Grand Tack evolution (see Fig. 
6) with the current distribution. Indeed, Minton and Mal- 
hotra (2011) showed that the current eccentricity distribu¬ 
tion can be achieved starting from a primordial distribution 
peaked at e ~0.3, similar to that produced by the Grand 
Tack evolution of Jupiter (Fig. 6).. They obtained this re¬ 
sult using secular resonance sweeping, but the basic result 
should hold also for a sudden enhancement of Jupiter’s ec¬ 
centricity. Nevertheless specific numerical simulations have 
never been done to demonstrate that the eccentricity distri¬ 
bution of asteroids at the end of the Grand Tack model can 
be transformed into the current distribution via the jumping- 
Jupiter evolution. The jump in Jupiter’s inclination due to 
its encounter(s) with another planet should also have par¬ 
tially reshuffled the asteroid inclination distribution, possi¬ 
bly reconciling the final distribution in Wetherill’s model 
(see Fig. 4) with the current one. However, the effects on 
the inclination during Jupiter’s jump seem less pronounced 
than those on the eccentricity (Morbidelli et al., 2010). 

The current inner edge of the asteroid belt is marked 
by the presence of the secular resonance between the pre¬ 
cession frequency of the perihelion of an asteroid and the 
fj(, planetary frequency. The resonance makes unstable all 
objects inside 2.1 at low to moderate inclinations, which 
tmncates the belt at its current edge. But before the im¬ 
pulsive separation between the orbits of Jupiter and Saturn, 
this resonance was much weaker (because the g 6 mode was 


less excited in the planetary system, the giant planets’ or¬ 
bits being more circular) and located away from the aster¬ 
oid belt. Thus, in principle the asteroid belt might have 
extended closer to Mars before the giant planet instability 
event. Bottke et al. (2012) showed that the destabilization 
of this extended belt -or ”E-belt”- could have dominated 
the formation of impact basins on the Moon, producing 12- 
15 basins over a 400 My interval. Given the age of the 
Orientale basin (usually estimated 3.7-3.8 Gy), this implies 
that the giant planet instability and the destabilization of 
the E-belt happened 4.1-4.2 Gy ago and was responsible 
for the production of the last dozen lunar basins known as 
the Nectarian and post-Nectarian basins. Earlier basins and 
craters would have to come from other sources, such as the 
planetesimals leftover from the terrestrail planet formation 
process. The existence of two populations of projectiles, 
namely the left-over planetesimals dominating the bom¬ 
bardment at early times and the E-belt asteroids dominat¬ 
ing the impact rate at a later epoch, should have produced a 
sawtooth-shaped bombardment history of the Moon (Mor¬ 
bidelli et al., 2012). The E-belt should also have caused 
a long, slowly decaying tail of Chicxulub-sized impacts 
on the Earth, possibly continuing until ~2 Gy ago. Evi¬ 
dence for this long tail in the time-distribution of impacts 
is provided by the existence of terrestrial impact spherule 
beds, which are globally-distributed ejecta layers created by 
the formation of Chicxulub-sized or larger craters (Johnson 
and Melosh, 2012): 10, 4, and 1 of these beds have been 
found on well-preserved, non-metamorphosed terrains be¬ 
tween 3.23-3.47-Ga, 2.49-2.63 Ga, and 1.7-2.1 Ga, respec¬ 
tively (Simonson and Glass, 2004; Bottke et al., 2012). 

Moreover, the escape to high-eccentricity orbits of bod¬ 
ies from the main belt and E-belt regions produced a spike 
in the impact velocities on main belt asteroids at the time 
of the giant planet instability. Thus, although the impact 
frequency on asteroids decreased with the depletion of 50% 
of the main belt population and 100% of the E-belt popu¬ 
lation, the production of impact melt on asteroids increased 
during this event because melt production is very sensitive 
to impact velocities (Marchi et al., 2013). For this reason, 
the impact ages of meteorites show a spike at 4.1 Gy like 
the lunar rocks, although for the latter this is due to a surge 
in the impact rate. 

A last consequence of the giant planet instability on the 
asteroid belt is the capture into its outer region of plan¬ 
etesimals from the trans-Neptunian disk (Levison et al., 
2009). Because Jupiter Trojans are captured in the same 
event, these last captured asteroids should have had the 
same source as the Trojans and therefore they should be 
mostly of taxonomic D- and P-types. The probability of 
capture in the asteroid belt is nevertheless small, so it is un¬ 
likely that an object as big as Ceres was trapped from the 
trans-Neptunian region in this event. 
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3.6. After the giant planet instability 

After the giant planet instability the orbits of all planets, 
giants and terrestrial, are similar to the current ones (within 
the range of semi major axis, eccentricity and inclination 
oscillation provided by the current dynamical evolution). 
Thus, the asteroid belt has finished evolving substantially 
under the effect of external events such as the giant planets 
migration or instability. 

The asteroid belt thus entered into the current era of evo¬ 
lution. Asteroids became depleted at the locations of un¬ 
stable resonances (mean motion and secular) on timescales 
that varied from resonance to resonance. In this way, the 
asteroid belt acquired its current final orbital structure. In 
this process, it is likely that another ~50% of the asteroids 
were removed from the belt, most of them during the first 
100 My after the giant planet instability (Minton and Mal- 
hotra, 2010). Combining this 50% with the 50% loss during 
the instability itself is the reason that we require that the pri¬ 
mordial depletion event (Wetherill’s model or Grand Tack) 
left a population of asteroids in the belt that was about 4x 
the current one. 

With the depletion of unstable resonances, the asteroid 
belt would have become an extraordinary boring place from 
the dynamical point of view. Fortunately, collisional break¬ 
up events keep refreshing the asteroid population, gener¬ 
ating dynamical families very rich in small objects, while 
non-gravitational forces, mostly the Yarkovsky effect (Bot- 
tke et ah, 2006) cause small asteroids to drift in semi major 
axis, eventually supplying new bodies to the unstable res¬ 
onances. This combination of collisional activity and non- 
gravitational forces allow the main asteroid belt to resupply 
and sustain in a quasi-steady state the intrinsically unstable 
population of Near Earth Objects. But this is the subject of 
another chapter. 

4. Conclusions and implications 

In this chapter, we have reviewed our current under¬ 
standing of the evolution of the asteroid belt, from a massive 
and dynamically quiet disk of planetesimals to its current 
state, which is so complex and rich from the points of view 
of both its dynamical and physical structures. 

According to this understanding, the asteroid population 
mainly evolved in two steps. There was an early event of 
strong dynamical excitation and asteroid removal, which 
left about 4 times the current asteroid population on orbits 
with a wide range of eccentricities and inclinations. This 
event may have been due to the self-stirring of a population 
of planetary embryos resident in the asteroid belt (Wether- 
ill, 1992) or to the migration of lupiter through the asteroid 
belt (the Grand Tack scenario; Walsh et ah, 2011). The sec¬ 
ond step occurred later, possibly as late as 4.1 Gy ago or 
~400 My after the removal of gas from the proto-planetary 
disk. At that time, the asteroid belt underwent a second 
dynamical excitation and depletion when the giant planets 
became temporarily unstable and their orbits evolved from 
an initial resonant and compact configuration to the current 


configuration. During this second event, the asteroid belt 
lost about 50% of its members. After this second event, the 
asteroid belt structure settled down with the progressive de¬ 
pletion at unstable resonances with the giant and terrestrial 
planets. Another 50% of the asteroid population was lost in 
this process, mostly during the subsequent 100 My. 

If the first evolutionary step was due to the Grand Tack 
migration of Jupiter, we expect that S-type asteroids formed 
more or less in situ (2-3 AU region); the C-type asteroids 
formed in the giant planet region (roughly 3-15 AU) and 
the P- and D-type asteroids formed beyond the initial loca¬ 
tion of Neptune (roughly 15-30 AU). The hot population of 
the Kuiper belt, the scattered disk and the Oort cloud would 
also derive from the same trans-Neptunian disk (Levison et 
ah, 2008; Brasser and Morbidelli, 2013). There is a growing 
consensus that the cold Kuiper belt f42—45 AU) is primor¬ 
dial and born in situ (Petit et ah, 2011; Parker et ah, 2011; 
Batygin et ah, 2011; Fraser et ah, 2014). Thus, the cold 
Kuiper belt objects should not have any correspondent in 
the asteroid belt. 

If instead the first step was due to the self-stirring of resi¬ 
dent embryos as in Wetherill’s model, we expect that S-type 
asteroids formed in the inner part of the belt, C-type aster¬ 
oids in the outer part and that no asteroids sample the plan- 
etesimal population in the giant planet region. The origin of 
P- and D-type asteroids would be the same as above. Thus, 
deciding which of these two models is preferable requires 
understanding better the nature of C-type asteroids, their 
water content, their similarities and differences with comets 
and of comets among themselves. This may not be an 
easy job to do. The population of main belt comets (aster¬ 
oids showing cometary activity such as 133P/Elst-Pizarro 
or 238P/Read) and their relationship with the parent bodies 
of carbonaceous chondrites is key in this respect. If it turns 
out that main belt comets are consistent with carbonaceous 
chondrites in terms of isotope composition (mostly for H, 
N and O) then this will argue that carbonaceous chondrites 
are just the rocky counterpart of bodies much richer in wa¬ 
ter/ice than meteorites themselves. This would imply that 
C-type asteroids formed beyond the snowline, thus presum¬ 
ably in the vicinity of the giant planets. If instead the main 
belt comets are not related to carbonaceous chondrites, but 
are more similar to comets from their isotope composition 
(it should be noticed that even though comet Hartley II has a 
chondritic water D/H ratio it has a non-chondritic 14 N/ 15 N 
ratio), then this would argue for their injection in the belt 
from the cometary disk and suggest that the parent bodies of 
carbonaceous chondrites are not so water-rich and therefore 
they formed somewhat closer to the Sun than the snowline. 

Whatever the preferred scenario for the first depletion 
and excitation of the asteroid belt, it is clear that the asteroid 
population must have suffered in the first hundreds of My as 
much collisional evolution as over the last 4 Gy. However, 
all asteroid families formed during the early times are not 
identifiable today because the dynamical excitation events 
dispersed them (and possibly depleted them) too severely. 
The presence of metallic asteroids not associated with a 
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family of objects of basaltic or dunitic nature, as well as the 
existence of rogue basaltic asteroids such as 1459 Magnya 
should therefore not be a surprise. The only families that are 
preserved are those that formed after the last giant planet in¬ 
stability event and that have not been made unrecognizable 
by subsequent collisional evolution and Yarkovsky drift; 
thus they are either relatively young or big. 

In this chapter, we have also examined several other as¬ 
teroid excitation and depletion scenarios, most of which 
have serious difficulties in reconciling their predicted out¬ 
comes with observations. We have done this not just for 
historical completeness, but also to illustrate the critical 
constraints on putative alternative scenarios of Solar Sys¬ 
tem evolution. For instance, numerous studies on the pos¬ 
sible in-situ formation of extrasolar super-Earths close to 
their host stars assume a large pile-up of drifting material 
of various sizes, from grains to small mass embryos, in the 
inner part of the prtoplanetary disks (Hansen, 2014; Chat- 
terjee and Tan; 2014; Boley et al„ 2014). By analogy, these 
models could be used to suggest that the outer edge of the 
planetesimal disk at 1 AU, required to form a small Mars, 
was due to the same phenomenon rather than to the Grand 
Tack migration of Jupiter. However, from what we reported 
in this chapter, we think that the asteroid belt rules out this 
possibility. In fact, the inward migration of small planetesi- 
mals (due to gas drag) and large embryos (due to disk tides) 
could explain the pile-up of solid mass inside 1 AU and the 
mass deficit of the asteroid belt, but not the asteroids’ orbital 
distribution. In the absence of the Grand Tack migration of 
Jupiter, we showed in section 3 that the only mechanism 
that could give to the belt an orbital structure similar to the 
observed one is Wetherill’s model of mutual scattering of 
resident embryos. But if this was the case, then the mass 
distribution could not be concentrated within 1 AU because 
a massive population of embryos is required in the main belt 
region. Thus, at the current state of knowledge (which may 
change in the future), only the Grand Tack scenario seems 
able to explain the required mass concentration to make a 
small Mars. 

In summary the asteroid belt remains the population of 
choice to test old, current and future models of solar system 
evolution. 
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